Oscillatory Flows Induced by Microorganisms Swimming in Two-dimensions 
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We present the first time-resolved measurements of the oscillatory velocity field induced by swim- 
ming unicellular microorganisms. Confinement of the green alga C. reinhardtii in stabilized thin 
liquid films allows simultaneous tracking of cells and tracer particles. The measured velocity field 
reveals complex time-dependent flow structures, and scales inversely with distance. The instanta- 
neous mechanical power generated by the cells is measured from the velocity fields and peaks at 15 
fW. The dissipation per cycle is more than four times what steady swimming would require. 



PACS numbers: 47.63. Gd, 47.15. G-, 47.61.-k, 87.16.Qp 

Understanding the hydrodynamics of swimming mi- 
croorganisms can give insight into complex biological pro- 
cesses such as predator/prey interactions and flagellar 
mechanics [JJ [2] ■ A variety of mathematical models have 
been proposed to describe interactions between swim- 
ming cells and their environment [3H5]- However, exper- 
imental studies of velocity fields are rare [5] , and oscilla- 
tory near-field hydrodynamics have not been investigated 
experimentally despite suggestion of their importance by 
some numerical studies [7j |8] . 

The flows induced by swimming organisms are respon- 
sible for biogenic mixing [S] , and microorganisms in par- 
ticular enhance the transport properties of suspensions 
[H [TU] . Active particles and microorganisms also affect 
the rheological properties of suspensions, where the ef- 
fective viscosity is increased in the case of "pullers" (e.g. 
Chlamydomonas reinhardtii) and decreased for "pushers" 
(e.g. E. coli) [TT1Q3. Hydrodynamic interactions leading 
to coherent motion have been studied for concentrated 
bacterial cells [15] . but not for algal cells. 

In many of these phenomena, the oscillatory swimming 
motion of organisms may be important. Furthermore, 
the concentration of swimmers may be sufficiently high 
that organisms are within several radii of one another, 
close enough that near-field effects should be considered. 
Understanding the resulting time-dependent flow fields 
is required for the interpretation of interactions between 
swimming cells and their environments. 

We present the first time-resolved measurements of the 
flow field around oscillatory swimming microorganisms 
(C. reinhardtii). This is accomplished through simul- 
taneous tracking of swimming cells and passive tracer 
particles within quasi-two-dimensional (quasi-2D) envi- 
ronments (» 15 /Ltm thick liquid films) |14j . where direct 
optical access and long observation times are achieved. 
These measurements reveal complex near-field flow struc- 
tures that evolve throughout the beat cycle, and provide 
a richer picture than time-averaged flow fields and simpli- 
fied models. The instantaneous mechanical power trans- 



ferred by flagella to the fluid is measured via the viscous 
dissipation. These observations carry important impli- 
cations for the interpretation and modeling of transport 
processes, locomotion, and flagellar mechanics. 

C. reinhardtii is a single-celled swimming alga with 
a 7-10 /im diameter body |15j . It typifies many biflag- 
ellated microorganisms with two 10-12 /im long flagella 
located at the anterior of the organisms. The flagella 
beat at 50-60 Hz propelling the cell body in an oscilla- 
tory manner (see video [TB]) with a mean velocity of 100- 
200 /im/s. The flagella are actuated along their length 
by dynein motors and take on asymmetric conformations 
during the power and recovery strokes to overcome the 
time reversibility of low Reynolds number flows. 

Wild-type C. reinhardtii (ccl690) is grown in minimal 
media (Ml) on a light cycle (16 hr bright/8 hr dark) to 
ensure uniform size distribution and motility (15] . An 
adjustable wire- frame device [12] is used to stretch a 2 /il 
droplet of the cell suspension into a square film with side 
length L = 6 mm and measured thickness h = 15 ± 2 /j,m. 
This thickness is sufficient to achieve 2D hydrodynamics 
for cells with radius R « 3.5 /jm [14) . A trace amount 
of gentle surfactant (Tween 20) stabilizes the film with- 
out harming the cells. Polystyrene microspheres (1 /nm 
diam., Thermo Scientific) are passivated with BSA to 
prevent adhesion to flagella [17] and incorporated with 
the cells prior to stretching the film. 

The device is mounted on an upright microscope where 
the suspensions are observed with a 40x 0.75 NA objec- 
tive under bright field illumination with red light (> 610 
nm) to prevent phototaxis [TS] ■ A high-speed digital cam- 
era captures the motion over about 3000 frames at 50 or 
500 fps. The cells and tracer particles are segmented 
from one another based on size and tracked in time using 
a predictive algorithm [TU] . 

The thin liquid film ensures that the algae are co- 
incident with the focal plane of the objective and pre- 
vents distortion from out-of-focus particles. This allows 
detailed observation of the cells (mean speed Uq = 134 
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FIG. 1: (color online) Probability density function (PDF) of 
flagellar beat frequencies, /, measured from unicellular alga 
(C. reinhardtit) swimming in quasi-2D liquid films (/ = 53 ±5 
Hz). Inset: velocity oscillations of a single swimming cell 
measured at 500 fps (red •), where the maximum velocity is 
four times the mean value (solid blue line, 134 fim/s). 



fim/s) for up to about 8 s at 50 fps. Using high-speed 
imaging (500 fps), oscillations of the cell body, U(t), be- 
come evident [Fig. [T] (inset)]. The peak forward velocity 
is four times the mean value and can be negative dur- 
ing the recovery stroke. The probability density function 
(PDF) of the beat frequency / for « 100 organisms is 
shown in Fig. [I] The oscillations have a fairly narrow 
distribution with a mean frequency / = 53 ± 5 Hz corre- 
sponding to a beat cycle period T = 1/f = 18.9 ms. 

By simultaneously tracking swimming cells and passive 
tracer particles (see video [Ej), we measure the velocity 
field induced by the cells by translating and rotating the 
instantaneous tracer particle measurements to a common 
coordinate system based on cellular orientation. The 
beat-cycle averaged velocity field is measured at 50 fps, 
and swimmer trajectory segments with large curvature 
or irregular beat frequencies are excluded [20]. The re- 
sulting time-averaged velocity field is shown in Fig. [2ja), 
where the swimmer motion is to the right and solid lines 
are instantaneous streamlines [2T]. The velocity field is 
calculated at 3 /iia resolution but represented at 6 /im 
resolution using log scaling to show the fine flow features. 

While the general shape of the flow field has some 
qualitative similarities to a force dipole ( "stresslet" ) ap- 
proximation [5], several unexpected features are appar- 
ent. The location of the hyperbolic stagnation point far 
away from the anterior of the cell is surprising, since it 
typically is thought to be located between the centers of 
drag (body) and thrust (flagella). Two strong vortices 
are visible lateral to the organism (associated with the 
flagella), while two weaker vortices appear far from the 
cell body beyond the separatrix of the hyperbolic point. 

In the quasi-2D liquid film, we expect a longer-range 
hydrodynamic disturbance compared to a 3D environ- 
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FIG. 2: (color online) (a) The time-averaged velocity field 
around a swimming C. reinhardtii (black disc) in the lab 
frame, where the direction of travel is to the right toward 
the hyperbolic stagnation point (green ♦). Solid (red) lines 
are instantaneous streamlines and velocity vectors are shown 
on a log scale [21] • (b) The fluid velocity magnitude in vari- 
ous directions away from the cell demonstrates the predicted 
u ~ r _1 scaling for a force dipole in 2D. Local minima corre- 
spond to stagnation points encountered in some directions. 



ment |14j . The magnitude of the fluid velocity normal- 
ized by the mean swimmer speed, u/Uq, where u = |u|, is 
shown in Fig. [2jb) as a function of the radial distance 
from the organism, r/R, in various directions. Toward 
the posterior, the fluid speed scales nearly as u ~ r -1 up 
to 20 cell radii away, consistent with a force dipole in 2D 
(compared to u ~ r~ 2 in 3D) [23] ■ The slight deviation 
from u ~ r _1 scaling is probably due to minor 3D effects. 
In the lateral direction, the velocity magnitude passes 
through a local minimum when traversing the parabolic 
stagnation points (vortex), before recovering u ~ r _1 scal- 
ing in the far field. Similarly, in front of the cell, the fluid 
speed decreases rapidly near the hyperbolic stagnation 
point 7-8 radii from the swimmer. 

The time-averaged velocity field in Fig. (2^a) shows the 




FIG. 3: (color online) A time sequence of the velocity field evolution throughout the beat cycle of C. reinhardtii (period T = 18.9 
ms) along with the position of the hyperbolic stagnation point (green ♦). Insets show cell speed and beat cycle phase (lower 
left, see Fig. [4jb) for details), and approximate flagellar shape (lower right) measured from high-speed video, (a) Early in the 
power stroke, the velocity field resembles a (negative) force dipole. (b) At the peak of the power stroke, the vortices lateral to 
the organism strengthen and sweep toward the posterior, (c-e) The vortices then shift to the anterior as the power stroke is 
completed, (f) At the peak of the recovery stroke, the flow field again takes the shape of a dipole, but with opposite sign. The 
recovery stroke velocity field (f) is weaker than the forward stroke, but is enhanced by the log scaling used in all panels [22 . 



limitations of some current swimmer models, and also 
raises several important questions. For example, how 
does the velocity field evolve throughout the oscillatory 
flagellar beat cycle? Using high-speed imaging (500 fps), 
we measure the instantaneous swimmer phase, and iden- 
tify tracer particles at corresponding times in the beat 
cycle. Velocity fields are constructed from tracer veloci- 
ties at each phase of the oscillation with resolution T/15. 

A time series of the velocity field during the beat cy- 
cle is shown in Fig. [3j where the cell is again swimming 
to the right (see video [16 ). Insets show the swimmer 
speed and phase (lower left) and the approximate flagel- 
lar position (lower right). At the beginning of the power 
stroke [Fig. [j^a)], the velocity field is neatly divided into 
four symmetric quadrants with the hyperbolic stagna- 
tion point located slightly forward from the body. This 
is in line with conventional force dipole swimmer mod- 
els [5]. As the fiagella move toward the posterior and 
the power stroke peaks, the vortices lateral to the organ- 
ism strengthen [Fig. |§fb)] and eventually shift across the 



body to the anterior side [Fig. |3fc-e)]. After the power 
stroke is completed and the recovery stroke begins, the 
fiagella extend out in front of the organism, and the cell 
velocity becomes negative. The flow shown in Fig. [3^f ) is 
qualitatively reversed from Fig. [3]ja), changing the sign 
of the dipole. The instantaneous flow field generated by 
an oscillatory swimmer such as C. reinhardtii is complex 
and highly time dependent. 

In Stokes flows, all of the mechanical energy generated 
by a swimmer for locomotion is rapidly dissipated by the 
fluid. The viscous dissipation per unit volume is given by 
$ = 2(T:r) where T is the rate of strain tensor. Thus, the 
total power transfered to the fluid at any instant by the 
organism's fiagella is P = J [i^hdA, where \i is the fluid 
viscosity and dA is a differential area element of the film. 
The phase-averaged velocity fields shown in Fig. [3] allow 
us to calculate the instantaneous mechanical power out- 
put, P osc (t), throughout the oscillatory beat cycle [Fig. 
HI a)]. The peak power output (« 15 fW) occurs during 
the power stroke and corresponds to the maximum in- 
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FIG. 4: (color online) The mechanical power dissipated by 
the organism is calculated from the velocity fields, (a) The 
instantaneous dissipation P OS c(i) is maximal when the swim- 
mer speed is highest. The dissipated power time-averaged 
over the beat cycle {P osc (t)) is more than four times the dis- 
sipation measured from the time-averaged velocity field (i.e. 
(Pose (t))/Pmcan flow ~ 4.2). (b) The measured power scales 
with the square of the cell body speed, P osc ~ U 2 . Inset: 
mean cell body velocity as a function of beat cycle phase. 



stantaneous speed of the cell body. A secondary local 
maximum also occurs at the peak speed of the recov- 
ery stroke [see Fig. |4|b) (inset)]. Due to the oscillatory 
swimming, the average mechanical power output over the 
beat cycle, (P osc (i)), is more than four times the power 
computed from the time-averaged velocity field in Fig. 

[fa) (i.e. (P osc (*))APmean flow ~ L t.^) . 

The mechanical power scales with the cell body speed 
[Fig. [4^b) (inset)] as P ~ U 2 with deviations at small U 
due to the finite accuracy of measuring tracer velocities 
[Fig. [4](b)] . This behavior might be expected for drag- 
based thrust at low Reynolds number, since the forces 
scale linearly with the body speed F ~ U, and the power 
output of the organism is P ~ FU ~ U 2 [H] . 

This work represents an important step toward un- 
derstanding the locomotion of swimming microorganisms 
and collective behaviors for which hydrodynamics are re- 
sponsible. Measurements of the time-averaged velocity 
field around C. reinhardtii in quasi-2D thin liquid films 
demonstrate the predicted u ~ scaling of the fluid ve- 



locity with detailed resolution of the near-field flow fea- 
tures. Also, high-speed imaging allows for phase-resolved 
measurements of the velocity field throughout the flagel- 
lar beat cycle, revealing complex underlying flow struc- 
tures that evolve in time. From these measurements, 
we demonstrate that the mechanical power dissipated by 
these swimming microorganisms is more than four times 
what the mean velocity field predicts. This has important 
consequences regarding the input of chemical energy re- 
quired to power oscillatory swimming of flagellated cells. 
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